Multiscale Modeling of the Active Site of [Fe] Hydrogenase: The H2 Binding Site in Open and Closed Protein Conformations

A series of QM/MM optimizations of the full protein of [Fe] hydrogenase were performed. The FeGP cofactor has been optimized in the water‐bound resting state ( 1 ), with a side‐on bound dihydrogen ( 2 ), or as a hydride intermediate ( 3 ). For inclusion of H₄ MPT in the closed structure, advanced multiscale modeling appears to be necessary, especially to obtain reliable distances between CH‐H₄MPT⁺ and the dihydrogen (H₂) or hydride (H^?) ligand in the FeGP cofactor. Inclusion of the full protein is further important for the relative energies of the two intermediates 2 and 3 . We finally find that hydride transfer from 3 has a significantly higher barrier than found in previous studies neglecting the full protein environment.     

QM/MM modeling of the hydroxylation of the androstenedione substrate catalyzed by cytochrome P450 aromatase (CYP19A1)

CYP19A1 aromatase is a member of the Cytochrome P450 family of hemeproteins, and is the enzyme responsible for the final step of the androgens conversion into the corresponding estrogens, via a three‐step oxidative process. For this reason, the inhibition of this enzyme plays an important role in the treatment of hormone‐dependent breast cancer. The first catalytic subcycle, corresponding to the hydroxilation of androstenedione, has been proposed to occur through a first hydrogen abstraction and a subsequent oxygen rebound step. In present work, we have studied the mechanism of the first catalytic subcycle by means of hybrid quantum mechanics/molecular mechanics methods. The inclusion of the protein flexibility has been achieved by means of Free Energy Perturbation techniques, giving rise to a free energy of activation for the hydrogen abstraction step of 13.5 kcal/mol. The subsequent oxygen rebound step, characterized by a small free energy barrier (1.5 kcal/mol), leads to the hydroxylated products through a highly exergonic reaction. In addition, an analysis of the primary deuterium kinetic isotopic effects, calculated for the hydrogen abstraction step, reveals values (～10) overpassing the semiclassical limit for the C? H, indicating the presence of a substantial tunnel effect. Finally, a decomposition analysis of the interaction energy for the substrate and cofactor in the active site is also discussed. According to our results, the role of the enzymatic environment consists of a transition state stabilization by means of dispersive and polarization effects. © 2015 Wiley Periodicals, Inc.     

Nucleic acid reactivity: Challenges for next‐generation semiempirical quantum models

Semiempirical quantum models are routinely used to study mechanisms of RNA catalysis and phosphoryl transfer reactions using combined quantum mechanical (QM)/molecular mechanical methods. Herein, we provide a broad assessment of the performance of existing semiempirical quantum models to describe nucleic acid structure and reactivity to quantify their limitations and guide the development of next‐generation quantum models with improved accuracy. Neglect of diatomic differential overlap and self‐consistent density‐functional tight‐binding semiempirical models are evaluated against high‐level QM benchmark calculations for seven biologically important datasets. The datasets include: proton affinities, polarizabilities, nucleobase dimer interactions, dimethyl phosphate anion, nucleoside sugar and glycosidic torsion conformations, and RNA phosphoryl transfer model reactions. As an additional baseline, comparisons are made with several commonly used density‐functional models, including M062X and B3LYP (in some cases with dispersion corrections). The results show that, among the semiempirical models examined, the AM1/d‐PhoT model is the most robust at predicting proton affinities. AM1/d‐PhoT and DFTB3‐3ob/OPhyd reproduce the MP2 potential energy surfaces of 6 associative RNA phosphoryl transfer model reactions reasonably well. Further, a recently developed linear‐scaling “modified divide‐and‐conquer” model exhibits the most accurate results for binding energies of both hydrogen bonded and stacked nucleobase dimers. The semiempirical models considered here are shown to underestimate the isotropic polarizabilities of neutral molecules by approximately 30%. The semiempirical models also fail to adequately describe torsion profiles for the dimethyl phosphate anion, the nucleoside sugar ring puckers, and the rotations about the nucleoside glycosidic bond. The modeling of pentavalent phosphorus, particularly with thio substitutions often used experimentally as mechanistic probes, was problematic for all of the models considered. Analysis of the strengths and weakness of the models suggests that the creation of robust next‐generation models should emphasize the improvement of relative conformational energies and barriers, and nonbonded interactions. © 2015 Wiley Periodicals, Inc.     

Mechanism of proteolysis in matrix metalloproteinase‐2 revealed by QM/MM modeling

The mechanism of enzymatic peptide hydrolysis in matrix metalloproteinase‐2 (MMP‐2) was studied at atomic resolution through quantum mechanics/molecular mechanics (QM/MM) simulations. An all‐atom three‐dimensional molecular model was constructed on the basis of a crystal structure from the Protein Data Bank (ID: 1QIB), and the oligopeptide Ace‐Gln‐Gly～Ile‐Ala‐Gly‐Nme was considered as the substrate. Two QM/MM software packages and several computational protocols were employed to calculate QM/MM energy profiles for a four‐step mechanism involving an initial nucleophilic attack followed by hydrogen bond rearrangement, proton transfer, and C? N bond cleavage. These QM/MM calculations consistently yield rather low overall barriers for the chemical steps, in the range of 5–10 kcal/mol, for diverse QM treatments (PBE0, B3LYP, and BB1K density functionals as well as local coupled cluster treatments) and two MM force fields (CHARMM and AMBER). It, thus, seems likely that product release is the rate‐limiting step in MMP‐2 catalysis. This is supported by an exploration of various release channels through QM/MM reaction path calculations and steered molecular dynamics simulations. © 2015 Wiley Periodicals, Inc.     

基础理论纵贯化学学科:吉林大学“统计力学与分子模拟”课程建设*

化学的基础理论的演进不断重塑着化学学科。传统的教学体系所依据的理论范式落后于现代化学理论的发展。吉林大学化学学院面向本科学生开设了“统计力学与分子模拟”课程。针对吉林大学本科学生的特点,精心安排了教学内容,通盘考虑知识体系的衔接,优化了教学方法。同时,淬炼了课程中的思政元素,从而在思想、道德、知识等3个方面进一步提升学生的综合素质。     

Error and stability analysis of numerical solution for the time fractional nonlinear Schrödinger equation on scattered data of general‐shaped domains

In present work, a kind of spectral meshless radial point interpolation (SMRPI) technique is applied to the time fractional nonlinear Schrödinger equation in regular and irregular domains. The applied approach is based on erudite combination of meshless methods and spectral collocation techniques. The point interpolation method with the help of radial basis functions is used to construct shape functions which play as basis functions in the frame of SMRPI. It is proved the scheme is unconditionally stable with respect to the time variable in and also convergent by the order of convergence , . In the current work, the thin plate spline are used as the basis functions and to eliminate the nonlinearity, a simple predictor‐corrector (P‐C) scheme is performed. It is shown that the SMRPI solution, as a complex function, is suitable one for the time fractional nonlinear Schrödinger equation. The results of numerical experiments are compared to analytical solutions to confirm the reliable treatment of these stable solutions. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1043–1069, 2017     

考虑底部隆起的浅埋隧道围岩压力计算分析

针对浅埋隧道开挖后底部隆起变形现象,应用极限分析上限法构造了考虑底部隆起变形的围岩压力计算模型,结合线性Mohr-Coulomb准则等推导出极限围岩压力的理论表达式。通过约束条件将围岩压力的计算转化为数学中的最优化问题,编制程序进行了优化计算。将计算结果与工程实测数据及文献计算结果进行了对比,验证了当前方法的可靠性。同时,指出在运用极限分析法处理浅埋隧道围岩压力问题过程中应将隧道底部一同考虑,对隧道底部支护会对整个围岩压力产生影响,有助于隧道结构的整体稳定。研究可以为浅埋隧道的开挖、支护提供一定的理论指导.     

Effects of Structural Variations on the Cellular Response and Mechanical Properties of Biocompatible,Biodegradable, and Porous Smectic Liquid Crystal Elastomers

The authors report on series of side‐chain smectic liquid crystal elastomer (LCE) cell scaffolds based on star block‐copolymers featuring 3‐arm, 4‐arm, and 6‐arm central nodes. A particular focus of these studies is placed on the mechanical properties of these LCEs and their impact on cell response. The introduction of diverse central nodes allows to alter and custom‐modify the mechanical properties of LCE scaffolds to values on the same order of magnitude of various tissues of interest. In addition, it is continued to vary the position of the LC pendant group. The central node and the position of cholesterol pendants in the backbone of ε‐CL blocks (alpha and gamma series) affect the mechanical properties as well as cell proliferation and particularly cell alignment. Cell directionality tests are presented demonstrating that several LCE scaffolds show cell attachment, proliferation, narrow orientational dispersion of cells, and highly anisotropic cell growth on the as‐synthesized LCE materials.